Equivalence between free and harmonically trapped quantum particles 
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It is shown that general solutions of the free-particle Schrodinger equation can be mapped onto 
solutions of the Schrodinger equation for the harmonic oscillator. This is done in such a way that 
the time evolution of a free particle subjected to a sudden transition to a harmonic potential can 
be described by a simple coordinate transformation applied at the transition time. This procedure 
is computationally more efficient than either state-projection or propagator techniques. A concate- 
nation of the map and its inverse allows us to map from one harmonic oscillator to another with a 
different spring constant. 



PACS numbers: 03.65.Ca, 03.65.Db, 03.65.Fd, 03.65.Gc 



I. INTRODUCTION 

It is known that laser beam modes in the parax- 
ial approximation are equivalent to the eigenmodes of 
the harmonic oscillator p], Q ■ The corresponding (one- 
dimensional, scalar) wave equation 



dx 2 



d 

9 (x, z) + ia —9 (x, z) = 



(1) 



arises as an approximation to Maxwell's equations |3|. In 
this optical context the parameter a = 4ir/\ is propor- 
tional to the (stationary, monochromatic beam's) wave 
number or inversely proportional to the light's wave- 
length A. 

The Schrodinger equation for a free particle in one spa- 
tial dimension x is given by 



d 

Q^2<j> (x, t) + ia —4> (x, t) = 



(2) 



where the parameter a = 2M/H contains the particle 
mass M and Planck's constant h, namely h = h/(2ir). It 
is obviously equivalent to eq. (U} above. 

Since the wave functions for free particles and those 
subjected to harmonic potentials factorizc with respect 
to their spatial coordinates we will only discuss the one- 
dimensional case. 

It has recently been realized that a simple coordi- 
nate transformation [2[ maps paraxial beams onto two- 
dimensional harmonic oscillator wave functions. Clearly, 
the same coordinate transformation can be used for the 
mapping of a free onto a harmonically trapped particle; in 
a more general setting this was noted before 0, @| • This 
case has added meaning since a non-adiabatic physical 
transition constitutes an experimentally implemcntablc 
sequence of changing environmental conditions. It turns 
out that modelling such a transition using the standard 
eigenstate-projection or quantum propagator techniques 
is more cumbersome. State projection typically leads to 



infinite sums over eigenfunctions which have to be trun- 
cated and are difficult to simplify, propagator techniques 
involve non-trivial integrals. The results reported here 
may therefore not only be of fundamental but also of 
technical interest. 

It is noteworthy that the discontinuities of the mapping 
of the time-coordinate in eq. ([3]) maps from a Schrodinger 
equation with a continuous to another with a discrete 
energy spectrum. 



II. MAPS 
A. Map from free to trapped case 

The coordinate transformations Q 



and t(r) = 6tan(rw) 



cos(tw) 

applied to the wave function mapping 

0(*(£,T),f(T)) 



where f(x,t;b) = 



/(*(£, r),t(r); 6) ' 



cxpl 4(fa+6a) 



(3) 



(4) 



(5) 
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(l + t 2 /6 2 )V4 ' 
yield the solution i/>(£,t) for the Schrodinger equation 

- |W r) + i* ^ (t r) + & r) = (6) 

of a harmonic oscillator with mass M, spring constant k 
and resonance frequency ui = y/k/M. 

Although this transformation maps t G [—00,00] onto 
r £ J— 7r/2, 7r/2), compare with the Gouy-phase of op- 
tics [21 , the periodicity arising through the use of trigono- 
metric functions meaningfully represents the oscillator's 
motion for all times r. 

In beam optics the confocal parameter b @ parame- 
terizes the strength of the beam's focussing and the cur- 
vature of its hyperbolic flow lines 0- Here, it serves as 



2 



a rescaling parameter of the transverse coordinate trans- 
formation and thus allows us to stretch or squeeze the 
width of the wave packet we want to map. 

In order to preserve the wave function normalization 
we have to determine the spatial coordinate stretching at 
the mapping time to = tq = 0. This yields the normal- 
ization factor M = (b w) 1 / 4 which the wave function has 
to be multiplied with. The specific choice 



&i = 1/cj = y/M/k 



(7) 



yields the natural mapping x — £ and Af = 1 which we 
will use from now on. 



B. Inverse map (trapped to free) 

The inverse to the coordinate transformations ([3]) are 



£ = 



and 



— arctan (uit) 



(8) 
(9) 



and go together with the wave function multiplica- 
tion <j) = ipf, inverting cq. 



C. Concatenated maps (trapped to trapped) 
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FIG. 1: Probability density P(x,t) ~ \(j>o(x,t;Q,p ) + 
4>o(x, t; 0, — po) | 2 of a free quantum particle in an equal su- 
perposition of two Gaussian states as described by eq. (|13|) ~). 
with h = 1, M = 1, <To = 3/2, xo = 0, and opposing mo- 
menta po = 4. The motion of the two half waves leads to 
interference at the origin. 
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Thc concatenated coordinate transformations from an 
initial harmonic trapping potential with spring con- 
stant k and wave function r; k), via the free particle- 
case 0, to a final harmonic potential with spring con- 
stant K and wave function r; k, K) is given by 



cos 2 (Or) + -g. sin 2 (fir) 



T = — arctan 



and *(^,T;fc,ii') 




tan (Or) , 



1&(3(£,r), T(r)) 
/(g(e,r),t(r);l) 
r),t(r);£)' 



(10) 
(11) 

(12) 



Here Q = K/M, and ^(^, r; fc, i^T) solves Schrodinger 
cq. ([6]) with A; substituted by K. 

As expected for a map from one harmonic oscillator 
to another, the inverse of the coordinate transforma- 
tions (fT0| and (fTTj) are given by the same functional ex- 
pressions with the quantities pertaining to one potential 
swapped with those of the other (k O K and u -*->■ f2). 



FIG. 2: Probability density P(x, i) of an initially free quan- 
tum particle which at time t = r = has the same state as the 
free particle displayed in Fig. [Tj and suddenly gets trapped by 
a potential with spring constant k — 5. This plot illustrates 
that the trapped particle is in a superposition state of two 
squeezed coherent states. 



III. AN EXAMPLE 

A well-known textbook example is the freely evolving 
Gaussian wave-packet with initial position spread ctq 



(f> (x,t;x ,po) 



1 



exp 



vq (xo ht — iMcr^x) 



tlOQ<j{t) 



x exp 
where a(t) = <tq + i 



(x — xq) i(Jo Mvo 2 t 



2<Jocr(t) 

th 



2ha(t) 



a M 



. (13) 
(14) 



Here Xq parameterizes spatial, and po = M v momentum 
displacement of the wave functions. If either of these two 
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quantities is non-zero the mapping onto a harmonically 
trapped state results in a state with oscillating center- 
of-mass. In general, although a wave function of Gaus- 
sian shape, this freely evolving wave packet will also not 
'fit' the width of the harmonic potential and therefore be 
squeezed. In short, the state of eq. (|13[) trapped in a har- 
monic potential becomes a squeezed coherent state Q. 
For a superposition of such states see Fig. Q] and for this 
superposition state being trapped consult Fig. [2] 

IV. CONCLUSION 

It is shown that general solutions of the free-particle 
Schrodinger equation can be mapped onto solutions of 
the Schrodinger equation for the harmonic oscillator us- 
ing a simple coordinate transformation in conjunction 
with a multiplication of the wave function by a suitable 



phase factor. This map is invertiblc and a concatenation 
of two such maps allows us to map from one harmonic os- 
cillator to another with a different spring constant. The 
simplicity of the approach described here makes it a tool 
of choice for the description of the wave function of a par- 
ticle experiencing instantaneous transitions from a free to 
a harmonically trapped state, the sudden release from a 
harmonic trap, or the sudden change of the strength of 
its harmonic trapping potential. 

The mapping introduced here is computationally more 
efficient than state-projection or propagator techniques 
and conceptually simpler than mapping techniques such 
as those used for supersymmetric potentials @, @ . It may 
even help with numerical calculations because it allows 
for the determination of the behaviour of wave functions 
trapped in a harmonic potential using the lower compu- 
tational overheads of wave functions in free space. 
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